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1. Introduction 


Regularization theory deals with the approximate solution of ill-posed operator 
equations 

nf) = 9 


in Hilbert or Banach spaces. In this paper we conhne ourselves to Hilbert spaces X and 
3^, and F maps from dom(F) G X to y. Let G dom(F) denote the exact solution 
and E y noisy data satisfying ||F(/1) — (?^|| < 6. One of the most prominent methods 
to obtain stable approximations to from such noisy data is Tikhonov regularization 


fS 
J a 


G argmin 

/gdom(F) 




( 1 ) 


A main question of regularization theory concerns the convergence \\fi — f^lx ^ ^ as 
5 —?■ 0 for appropriate choices of a = a{6,g^), and the rate of this convergence. To 
obtain such rates, additional assumptions on /I are required (see [12, Prop. 3.11]). For 
a long time such conditions were formulated as spectral source conditions of the form 
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yt = (^(F'[/^]*F'[/^])w (see [38] and numerous references therein) for some w & X and 
an index function ip, that is ip : [0, cx)) —)■ [0, cx)) is a continuous increasing function 
satisfying y9(0) = 0. Starting with [23] source conditions have been formulated in the 
form of variational inequalities 


V/ G dom(F) 


2 


f-f 


2 1 
< - 
X - 2 




( 2 ) 


with some parameter fd G (0,1], and this type of source conditions has become more 
and more popular in regularization theory (see e.g. [3-5,7,9,11,15-17,19,21,22,24-26, 
36,37,45]) due to the following advantages over spectral source conditions: 


Proofs based on variational source conditions tend to be much simpler than proofs 
based on spectral source conditions. E.g., for Tikhonov regularization and concave 
index function ip a simple argument by Grasmair [20] (see also [45, Thm. 3.3]) yields 
the convergence rate 

/3 




( 3 ) 


for an optimal choice of the regularization parameter a fulhlling —1/(2 q;) 
where d{—'ip) denotes the subdifferential of —ip. 


• For bounded linear operators in Hilbert spaces variational source conditions are not 
only sufficient, but even necessary for certain rates of convergence for Tikhonov 
regularization and other regularization methods in many cases (see [18]). For 
spectral source conditions this is only true in a supremum over /f in certain 
smoothness classes, but not for individual /h The result in [18] were derived for 
so-called approximate source conditions, which have been shown to be equivalent 
to (2) in [14,15]. 

• Variational source conditions do not involve the Frechet derivative F' of the forward 
operator. This reduces the regularity assumptions on F, but more importantly, it 
avoids restrictive assumptions on the relation of F and F' such as the tangential 
cone condition, which cannot be verihed for most interesting nonlinear applications, 
e.g. inverse scattering problems. 


• Variational source conditions can be used not only in a Hilbert, but also in a Banach 
space setting and for more general noise models and data hdelity terms (see e.g. 
[11,26,45]). 


However, so far variational source conditions could be verihed only in rather few cases: 
One option is to derive them from spectral source conditions, but then the variational 
approach does not yield additional information. For linear operators F and F penalties 
with respect to certain bases in the range of F*, characterizations of variational source 
condition have been derived in [4,7]. Moreover, reformulations of (2) with ip{x) = \/x 
for a phase retrieval and an option pricing problem were derived in [23]. The purpose of 
this paper is to show for a classical inverse scattering problem that a variational source 
conditions holds true under Sobolev-type smoothness assumptions. 
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The forward problem we consider is as follows: Given a refractive index n = 1 — f 
and one or several incident wave(s) u'" solving the Helmholtz equation = 0, 

determine the total £eld(s) u = u' + u® such that 

An + K^nu = 0 in (4) 

du"" . ^ f l\ 

-- iKU = u [ — ] as r = X —)• cx). (5) 

or \r^ / 

We will study inverse problems to recover the refractive index given measurements of 
scattered helds as formulated precisely in Section 2. Related problems occur in many 
applications in nondestructive testing, geophysical exploration, and x-ray imaging. 

The main tool of our analysis are geometrical optics solutions introduced in 
Section 3. The use of such functions is well established for deriving uniqueness and 
conditional stability estimates, which for convenience we write in the form 

fll/i-Mlt<V-(l|f(/i)-F(/2)llb (6) 

for certain (typically compact) smoothness classes K C dom(F), e.g. Sobolev balls. For 
the acoustic inverse medium scattering problem the hrst such estimate was established 
by Stefanov [42] using a very strong norm in the image space y and a logarithmic 
function of the form 

ip{t) = G(lnf“^)“^^ for some /i > 0. 

This estimate was improved in [29] by choosing y as an L? space and making the 
exponent fi explicit with p < 1. Recently, improved estimates were established in 
[31] where /i —?■ cx) as the smoothness exponent of the Sobolev ball K tends to oo. 
Lower bounds in [30, 40] show that such upper bounds are essentially optimal. The 
dependence on the wave number n was made explicit in [32, 34] leading to so-called 
Holder-logarithmic stability estimates with signihcantly improved stability for large k. 

To make this hrst paper on the verihcation of variational source conditions accessible 
to a larger audience and to keep the level of technicality as low as possible, we did not try 
to incorporate all of these recent improvements for stability estimates into our analysis. 

Let us compare variational source conditions (2) and conditional stability estimates 
(6) on the abstract level of operator equations. Obviously, (2) for all G iF implies 
(6) since we may assume w.l.o.g. that ||/i|| > II/ 2 II and choose f = fi and /I = / 2 . 
However, the reverse implication is not obvious since the case ||/|| < ||/1|| cannot be 
excluded and since (2) is required for all / in the larger set dom(F), not only f E K. 

It is interesting to have not only a stability estimate (6), but also a variational source 
condition (2) since the latter yields error bounds for reconstructions from noisy data 
obtained by Tikhonov regularization and other commonly used regularization methods. 
From conditional stability estimates one can also derive error bounds for the method 
of quasisolutions G argminjg^||F(/) — g^\\y (see [35]). However, this method is 
often difficult to implement and rarely used in practice (see e.g. [2] for a theoretical 
discussion for an inverse obstacle scattering problem). Moreover, the set K must be 
known explicitly, a typically unrealistic assumption, whereas Tikhonov regularization 
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with a-posteriori selection rules for the regularization parameter a can attain optimal 
rates of convergence over a large range of smoothness classes without a-priori knowledge 
of the smoothness of the solution. 

The plan for the remainder of this paper is as follows: The next section contains 
a precise formulation of our main results. Section 3 describes properties of geometrical 
optics solutions, and Sections 4 and 5 contain the proofs of our main results. 


2. Main results 


We assume that the contrast / = 1 — n is supported in the ball 53 (vr) where 53 (i?) : = 
{x G 

kl — R} for R > 0. Due to the physical constraints 3?(n) > 0 and 9^(n) > 0, 
the contrast / belongs to the set 

©:={/£ L“(E3): SJ(/) < 0, »(/) < 1, supp(/) C »(!r)} . 

It is well known that for all / G D a solution u to (4) exists, and it is unique due to (5). 
Corresponding inverse problem consist in recovering / from measurements of certain 
total helds u for different incident helds u\ We will discuss two such problems below. 
The hrst one is to reconstruct / from near field data 

Wf{x,y) = u^ix) +ul{x), {x,y) G d^{R) x d^{R), 

for some R> tt. That is we measure for each incident point source wave of the form 

^iK\x-y\ 

centered at y G d^{R) the corresponding total field on the sphere d^{R). Our aim is 
to solve the equation Ffif) = w for given data w, where 

F^-.V ^ L\d^{Rf), f^Wf 

is the near field operator that maps the contrast / to the Green’s function Wf of the 
pde. 

As preimage space X we choose the Sobolev space HlfiifiBij)) equipped with the 
norm 

1/2 


- E (i + hi'T \R'> 


where /(q) = {2Tr)~^^'^ f{x)e~^'^'^dx are the Fourier coefficients of / in the cube 

(t(7r) = (—TT, tt)^ with / extended by 0 outside of 53(7r). Our main result for near held 
data reads as follows: 

Theorem 2.1 (Variational source condition for near held data). Let | < m < s, 
s 7 ^ 2m -1-3/2 and tt < R. Suppose that the true contrast /I satisfies /I G D with 
ll/^llm < Cg for some Cg > 0. Then a variational source condition (2) holds true for 
the operator F^ with dom(Fn) := V D X, y = , [3 = 1/2, and given by 

—2fi 


finit) := A (ln(3 -ht ^)) 


m 


y := mm 


’ m + 3/2^ 
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where the constant A> depends only on m, s, Cg, k, and R. 

Actually, this theorem holds true for arbitrary /3 G (0,1), but the constant A 
depends on the choice of /3 (see Remark 4.2). For a discussion of the case s = 2m + 3/2 
see Remark 4.4. From the result cited in (3) we obtain a rigorous error bound for 
Tikhonov regularization without further assumptions on Fk: 


Corollary 2.2. Under the assumptions of Theorem 2.1 the error bound 




-V 


in terms of the noise level S holds true for the regularization scheme (1) if a is chosen 
such that ^ 

We could formulate further corollaries from regularization theory, in particular on 
a-posteriori choice of the regularization parameter. Moreover, we obtain the following 
stability estimate: 


Corollary 2.3. Let | < m < s, s 7 ^ 2m + 3/2 and n < R. Let fi and f 2 satisfy /j G "D 
with ||/j||r/s < Cg for i = 1,2 and some Cg > 0. Then the stability estimate 

ll/i - / 2 ll„» < 2\/I(hi (3 + ||r„(/i) - .,))■" (7) 

holds true with p and A as in Theorem 2.1. 


(7) differs from similar results in the literature discussed in the introduction in the 
use of a Sobolev norm rather than an or an norm on the left hand side. However, 
this could easily be accommodated for in proofs of previous stability results, so the 
relevance of this work is rather expressed in Corollary 2.2 than in Corollary 2.3. 


For the formulation of the second inverse problem recall that every solution u of 
the Helmholtz equation (4) fulhlling the Sommerfeld radiation condition (5) has the 
asymptotic behavior 

u{x) = +(a;) + - — + O > r = |x| —)■ cx) 

uniformly for all directions x = x/r G & := {x G |a:| = 1}, and u°° is called the far 
held pattern (see [10,39]). Often far held patterns u°°{-,d) corresponding to incident 
plane waves u]j^^{x) = propagating in direction d are considered as data. This 

leads to an inverse problem for the forward operator 

FfiPn A ^+2(6^ X ©2), f^u°°. 

Theorem 2.4 (Variational source condition for far held data). Let the assumptions 
of Theorem 2.1 be satisfied. Then the operator F} with dom(Ff) -.= T> A X and 
Y = L\G'^ X ©^) fulfills for all 0 < 9 < 1 a variational source condition (2) with 
j3 = 1/2 and given by 

V'f(t) := R (ln(3+ f^)) ^ 

with pi given as in Theorem 2.1, and a constant B > 0 depending only on m, s, Cg, k, 9, 
and R. 
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Once again we obtain a convergence rate via (3) and a stability estimate as corollary: 
Corollary 2.5. Suppose the assumptions of Theorem 2.4 hold true. 

(i) Then the minimizers of the Tikhonov functional (1) satisfy the error bound 


ft - /' 




Hr, 


< 


4/B (ln(3 + 5-2; 


—fiO 


t=46^ • 


(a) For all /i ,/2 G V satisfying ||/j||h'> < Cg for j = 1,2, s ^ 2m + 3/2 and some 
Cg > 0 the following stability estimate holds true: 

ll/i - Alii,™ < 2x/B (in (3 + ||f,(/i) - ff(/2)|IZ.Vx6.)))''‘'' ■ 


3. Bounding Fourier coefficients using geometrical optics solutions 


As mentioned in the introduction, the key tool of our proof are geometrical optics 
solutions. These functions are solutions of the perturbed Helmholtz equation (4) of 
the form 

M(a;) = e*^^(l+ n(x,C)), (8) 

where C G \ satisfies C • C = (where ( ■ f '.= for ^ C^). Such 

exponentially growing solutions were introduced by Faddeev in [13] and have been used 
to prove uniqueness of electrical impedance tomography and inverse scattering problems 
[8,41,43] as well as for stability estimates [1,29,31,32,34,42] for inverse scattering for 
space dimension d > 3. We recommend the textbook [39] or the monograph [10] for 
concise and self-contained introductions of the version of geometrical optics solutions 
used below and refer to the monograph [33] and the review [44] for numerous extensions 
and further references. 

A reason for the interest in geometrical optics solutions is the following lemma: 


Lemma 3.1 ([29, Lemma 3.2]). Let tt < R < R'. Then there exists a positive constant 
Cl (depending on k, R, and R') such that for all contrasts /i, /2 G H with corresponding 
near fields Wi and W 2 , and for all solutions Uj G H‘^(f3{R')) to Auj -|- = ri^fjUj the 

following estimate holds true: 


/ (/i - / 2 ) uiU 2 dx 

-^Q3(7r) 


< Cl llwi 


II'“i|Il2(33(K')) II'“2||l2(3s(/j/)) . 


Note that if we could choose Ui and U 2 as geometrical optics solutions with (/i = ^'2 
and Vi = V 2 = this lemma would immediately yield bounds on the Fourier coefficients 
of fi — f 2 - Even though such a choice is impossible for /i, /2 ^ 0, we will derive bounds 
on the small Fourier coefficients of fi — f 2 in Lemma 3.3 using bounds on v detailed 
below. 


Proposition 3.2 (Existence and norm estimate of geometrical optics solutions). Let 
R' > 71 and / G H. Then for all G with C • C = such that t := |Q'(C)| fulfills 
t > 2^2:^ II/|Iloo(r 3 ) there exists a function v{x, () such that u{x, () := (1 -1- v{x, ()) 
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belongs to solves the eguation Au + n?u = in 53(7?'), and satisfies the 

estimates 


|n(', C)llL2(g5('^/)) < ^ |l/|lLtx>(R3) 

(9) 

n(', C)IIl2(<8(_R')) ^ O2 

(10) 

«(gC)IIl2(^(k,)) < C2e^'l®(^)l 

(11) 


with a positive constant C 2 depending on k,R, and R'. 

Proof. For the mapping properties see [27, Lemma 5], the hrst estimate can be found 
in the proof of [28, Lemma 2.9]. To prove ( 10 ), insert the lower bound 2 ^^^ II/IIl°°(r3) 
for t into (9). As 

'e<-(l + w(a;,C)) 


K,C)II 


L2(3S(R')) 


< 






one obtains ( 11 ) by inserting ( 10 ). □ 

In the following the constant Mem is given by the Sobolev embedding theorem such 

that 

||/IL.(^)<Mem||/||H^ (12) 

for all / G i 7 ™(f 8 ( 7 r)). The following estimate is similar to the weaker result [29, Lemma 
3.3] where the factor ||/i — f 2 \\H^ on the right hand side is missing. A similar result 
can be found in [32, Lemma 2.1], but its proof relies on more sophisticated results on 
geometrical optics solutions. 

Lemma 3.3. Let Cm > 0, m > 3/2, tt < R < R' and fi and f 2 be contrasts with fj E R 
and II/jIIh"' < Cm with corresponding near field data Wj for j = 1,2. Define 


7?' 

fn := 2k^ - MemCm 


TT 


(13) 


with Mem defined as in (12). Let t > to and 1 < g < 2\/n? + R. Then there exists a 
constant C 3 > 0 depending only on m, R, 7?', and k such that for all 7 G satisfying 
I 7 I < g we have 


(/i - / 2 )( 7 )| < llwi - W2|lL2(a>8(/?)b + y ll/i - /2II 


Hm ■ 


Proof. For hxed 7 G choose two unit vectors di and d 2 in such that 7 • di = 
7 • (72 = (7i • (72 = 0. For t > R dehne 


Ct ■= -|7 + itdi + y ^2 + 


d 2 


Vt •= “2^ “ “ V ^ ~ ^ 

Then Ct,Vt ^ satisfy 

Ct + vt = -i, |^(Ct)l = 1^(701 > ^0, CfCt = hfVt = 
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Hence by Theorem 3.2 there exist geometrical optical solutions of the form 

ui{,x, Ct) = (1 + vi{x, Ct)), 

Mx, Vt) = (1 + V2ix, r]t)) , 

where Uj solves the equation Auj + = rc^fjUj in for j = 1, 2. It follows that 


(/i - / 2 ) ( 7 ) 


(27r)3/2 

1 


/23(7r) 


(/i-/2)(x)e-*"-"dx 


' {fl - f2){x)Ui{xXt)u2{x,T]t)dx - 

23(7r) J^{7v) ^ 


(27r)3/2 

(ni(a:,Ct) + V2ix, T]t) + Vi{xXt) V2ix,r]t)) 


{fi- f2){x) 


dx 


(14) 


The first integral on the right hand side of (14) can be bounded by Lemma 3.1 and 
by the norm estimate (11) for geometrical optics solutions for t > to: 


'2S(7r) 


(/i - 12 ){x)ui{xXt)u2{x,r]t)dx 


< C 3 \\Wi-W 2\\L2^g<s(Rf)e^^'^- 


Using (9) and (10) the second integral on the right hand side of (14) can be 
estimated by 




e-*^'"(/i 

- J 




< 

1/1- 

h\\L‘2{<3{7 

0 ) ( 


+ 11/1 

~ / 2 |Ilcx> 

(®(r 

<fllA 

~ /2|Lffm 



(16) 

□ 


Plugging (15) and (16) into (14) yields the assertion. 

To avoid another free parameter, we will set R' := 2R in the rest of this paper. 

4. Proof of Theorem 2.1 


We will prove the following equivalent formulation of the variational source condition 

(2): 


S (/',/*-/ 




S E ( 

V 76 Z 3 


1 + - /)(7) 


< 


1-/3 


(17) 


/^-/|L^+V^( FXf)-FXP) 


y 


This form has also been used in regularization theory (see e.g. [23]), but we preferred 
to formulate our main theorems using the form (2) since its relation to stability results 
is more obvious and since the form (2) is used in the proof of the convergence rate (3). 
We bound the high Fourier coefficients on the left hand side of (17) as follows: 
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Lemma 4.1. //^ > 0, 0 < m < s, p ^ Hp^{Ti)) and f G then 

i S 


3f? 


\7GZ3\Q3(e) 


1 + 


PiP {P - f)P) 


< 


P-f 


+ 2 




Q' 


2{ra—s) 


Proof. By the Cauchy-Schwarz inequality and Young’s inequality xy < 2x^ + for 
x, 1 / > 0 we have 


^1 E (1 + 


PP){P - f)P) 


< 

E 

( 1 +l 7 lT 

{p - f) ( 7 ) 

Z3\a3(e) 


1 

< - 

- 8 

E 

(i l 7 lT 

{p - /) ( 7 ) 




PP) ( 18 ) 


The hrst sum on the right hand side is bounded by |||/^ — To bound the other 

sum we use that /I is smoother than / to obtain 


< 


E {^+bfr\fKP - . ^ 

Z3\!B(e) 13- W 

Inserting these estimates into (18) completes the proof. 


E (i + hlTk’Wrse^'’”->||/' 




□ 


Obviously this bound is only useful with a proper choice of g as detailed below. To 
bound the low Fourier coefficients we use Lemma 3.3: 

Proof of Theorem 2.1. Given / G fl P we distinguish two cases: 

Case 1: ||/^ — /||irm > ACs- By the Cauchy-Schwarz inequality we have 




< 


P P -f 






P-f 


< - 

— A 


P-f 


(19) 




which clearly implies (17). 

Case 2: ||/1 — fWu”^ < ‘^Cs- Then < 5Cs, and hence we can apply Lemma 

3.3 with Cm = SC's and to as dehned in (13). Moreover, we choose R! = 27?, t > R, and 
1 < i? < 2p+ P and set 5 := \\wf — ujp\\L 2 (gfs{Rf) 1° obtain 


s( E (i-H7lT7h)(7-/)(7)') <C3 E |(i + l7lT/'W 

Vz3n'B(p) / 7ez3n58(e) 

f-P ) 


X 


X e'‘“S + - 


( 20 ) 
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The sum on the right hand side of (20) can bounded by the Cauchy-Schwarz inequality 
and Lemma 4.3 below with A = 2m — s as follows 


E 

7ez3n®(e) 


(i + i7iThM 


< 


\ -yei? 
< C^CsQ 


E (i + biT k'b) 


\ 7ez3n58(e) 


(i + IaI' 


2m—s 


( 21 ) 


with r = max{2m + 3/2 — s, 0}. Putting Lemma 4.1, (20), and (21) together yields 




< 


c^c^CsQ" 


/t - / + > + €3040.0““+ 

n t 

1 1 ^2r+2(s-m) ^ ^ 


f-f 


H-n 


( 22 ) 


< x + t; 


8 et"^ 


- f\\^^ + C3C4C',e^«V"5 + 2^7/ (l + eclcl) 


for all e: > 0 by Young’s inequality. Now we have to choose the free parameters t, g and 
e such that the right hand side of the last inequality is approximately minimal for given 
5 > 0 and the constraints in Lemma 3.3 are satisfied. We pick t, g and e such that 


12Rt = ln(3 + 6-^) = e := {12R)^. 


(23) 


As T+s—m > m+3/2 > 3, there exists t > 0 such that 2t > for all t >t. 

Let us strengthen the constraint t > to in Lemma 3.3 to t > to with to := max{to,t}. 
Then the constraint 2 \/+ R > 2t> (12i?t)^/’^'’"+^“”’'i = g in Lemma 3.3 holds true for 
all t > to, and p > 1 is satisfied as well. However, with the choice (23) the constraint 
t > to itself is only satisfied for 5 < ^max with ^max := (exp(12i?to) —3)“^/^ (or (5max := oo 
if exp(12i?to) < 3). The case 5 > (5max will be treated at the end of this proof. Plugging 
(23) into (22) yields 




1 

< - 
- 4 


< 


P - f + 2Cs (1 + ^^ 304 ) (ln(3 + 6 ^)) 
C 3 C 4 C ,(3 + (ln(3 + 


_27\-2m 


■f’ - f\\„~ + ^ (‘"(3 + 3-t) 


2 s — m, 
r+s —m 


for S < (fjnax with some constant A depending only on m, s, Cg, n, and R since the term 
in the second line tends to 0 faster than the last term in the first line as 5 \ 0. This 
shows (17) for 6 < (5 max . 

It remains to treat the case 6 > (5 max - By the Cauchy-Schwarz and Young’s 
inequality we have 




1 

< - 
- 4 

/*-/ 

2 

JJm 

P 

2 

Hm 

1 

< - 
- 4 

f-f 

AC's 

}jm * 



Hence, Theorem 2.1 holds true with A max{4, C^(ln(3 + d,„^))^'‘}. 


(24) 

(25) 
□ 
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Remark 4.2. Note that in the previous proof both contributions to the term 
^^11/^ ~/111/'" ill (li") originate from an application of Young’s inequality xy < 
^x^ + ||/^. Hence, by a proper choice of £ > 0 Theorem 2.1 can be shown for arbitrary 

/9 e (0,1). 

It remains to show the following lemma needed in the previous proof: 


Lemma 4.3. For all A ^ —3/2 there exists C 4 > 0 depending on A such that 

^ 7GZ3nQ3(e) 

for all p > 1 where r = max{A + 3/2, 0}. 


(26) 


Proof. For each 7 G we introduce a cube := {x G : ja: — 7|cx) Y 1/2} around 7 
with side length 1. Then U^gz ^7 = two different cubes intersect at most on a 

set of measure 0. Consider hrst the case that A > 0. As Y 53(p + \/3), it 

follows that n iB(p)} < vol + \/3)), so we can estimate 

E (l + lbl^ < + forp>l. 

7ez3n'B(p) 


Now consider the case A < 0. For 7 7 ^ 0 note that C.y fl fB(| 7 |) contains the ball with 
radius 1/4 centered at 7 — A sgn( 7 ), the volume of which is Y = 7r/48. Hence we obtain 

. A 


E (1 + 

'y&’Z^n<3{e) 


<1 + / E / (l + klAd^ 


< 


1 + — [' rUl + r^Y dr <1 + 192 [ 
V Jl/2 ^ A 

1 + 0 > A > -3/2, 

^ (-2A-3)23+2^ ’ ^ —3/2. 

This shows the lemma for all values of A 7 ^ —3/2. 


„2+2A 


1/2 


dr 


< 


□ 


Remark 4.4. It is easy to see that for A = —3/2 the left hand side of (26) grows 
logarithmically in A. Using this fact, it is straightforward to treat the case s = 2m+ 3/2 
in Theorem 2.1 in a similar way, which leads to a function //^n involving an additional 
logarithmic term. Since we believe that this upper bound is not optimal anyways, we 
did not consider it worth to include the computations for this special case. Note that if 
/I G for s = 2m + 3/2 then also /I G H^' for all s > s', and therefore Theorem 2.1 
still holds true with s replaced by s' for all s > s' > m. 


5. Proof of Theorem 2.4 

In this section we prove the extension of the result for near held data to far held data. 
The main tool is the following lemma relating near and far held data: 
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Lemma 5.1 ([29, Section 4]). Let R > ir, m > 3/2, Cm > 0 and 0 < 9 < 1. Then 
there exist constants cv, g, 5max > 0 such that for any two contrasts fi, f 2 ^T) satisfying 
< Cm we have 


lk2-M^i|li2(a^(2K)2) < ^"exp - 


— In 


Un 


u 


1 IIl2(62x62) 


ug 


(27) 


< ^max where Wj and denote near and far field scattering data 

for fj, J = 1,2. 


Proof of Theorem 2.4- We distinguish the same cases as in the proof of Theorem 2.1: 
Case 1: ||/^ — > ACg- This case can be treated as in the proof of Theorem 2.1. 

Case 2: < dC*. Then < SC*. Hence Lemma 5.1 can be 

applied with Cm = SC's. Setting 5 := ||Cf(/l) - Tl(/)||L 2 (e 2 x 62 ) and : = 

exp(—(—ln(\/t) + ln(a;p))^), it follows from Theorem 2.1 and the monotonicity of 
that the variational source condition (2) holds true with fjit) = ' 0 n(<^(t)) if ^ < ^max- 
Bounding '0n(t) < A{\\it~^)~‘^r for f < 1 we obtain 

' 0 n(<y( 5 (t)) < A (-ln(\/t) + In(pcu)) -Inp^^ for \/t < min |(5max, ^| • 

Hence, it is easy to see that there are constants H > 0 and hmax G (0, min/h max . |}] such 
that 

'0n(95(i)) < H(ln(3 + for \/f < ^max- 

This shows (2) for 5 < ^max- The case 8 > ^max can be treated as in the proof of Theorem 
2.1, see eq. (24). □ 


6. Conclusions and Outlook 

We have established a source condition in the form of a variational inequality for three- 
dimensional acoustic medium scattering problems which implies logarithmic rates of 
convergence of Tikhonov regularization under Sobolev smoothness assumptions as the 
noise level tends to zero. Our proof is based on geometrical optics solutions, a technique 
that has been developed to show uniqueness and stability results for a number of 
parameter identification problems in differential equations. 

Our results may be extended in several directions: One may study other differential 
equations, in particular other inverse scattering problem such as electromagnetic 
or elastic scattering problems ([28]), electrical impedance tomography ([!]) or two- 
dimensional acoustic scattering problems ([ 6 ]). Moreover, analogs to improved stability 
estimates with explicit dependence on the wave-number and exponents tending to cx) 
with the Sobolev smoothness index would be of interest ([32]). Finally, variational 
source conditions based on Banach norms would be desirable. So far we do not see 
a general scheme for establishing variational source conditions based on conditional 
stability estimates, it rather seems that each extension requires new ideas. 
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Even beyond the technique of geometrical optics solutions it seems a worthwhile 
endeavor to study which of the large variety of conditional stability estimates for 
various forward operators and smoothness classes can be sharpened to variational source 
conditions. 
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